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1.0  INTRODUCTION 

Among  the  many  effects  to  consider  in  modeling  environmental  acoustics  is  the 
phenomenon  of  forward  scattering.  This  phenomenon  is  often  referred  to  as  out-of-plane 
scattering.  For  a  rough  surface,  it  will  result  in  a  distribution  of  energy  over  some  angular 
region,  as  well  as  a  loss  of  coherence. 

The  purpose  of  this  document  is  to  develop  theories  that  will  address  various  issues 
pertinent  to  the  problem.  A  theoretic  outline  drawn  from  various  sources  is  presented  in  a 
consistent  fashion.  The  citation  of  these  sources  is  in  itself  useful.  In  addition,  some 
details  of  the  numeric  implementation  are  discussed. 


2.0  THEORETICAL  BACKGROUND 

Consider  the  problem 


V2px  +  k2pr  =  -/( r) 


(1) 


for  the  pressure  wave  Pt  with  wavenumber  k  from  a  distributed  source  /.  The  solution  in 
terms  of  the  Green’s  function  G(rjR),  where 


V2G  (r|R)  +  k2G  (r|R)  =  -6  (r  -  R)  ,  (2) 

is 


Pt(  r) 


/  j  J  /(R)G'(r|R)d3R  + 

J  J  |c(r|R)  ~pt(R)  —  Pt(R)  ~G(rjR) 


d2R. 


(3) 


For  waves  radiated  into  a  unbounded  medium. 


G(r|R) 


1  ei*|r-R| 
47r  |r  —  Rj 


(4) 


In  the  scattering  problem,  the  volume  integral  is  zero  and  the  surface  integral  is  over  the 
sphere  at  infinity  and  the  object  of  interest.  The  integral  over  the  sphere  at  infinity  pro¬ 
duces  the  incident  wave  po .  Thus 


pr( r)=Po(r)  +  J  ^[pr(F  )^G(r|R)-G(r|R)^pr(R) 


dS  , 


(5) 


which  is  an  integral  equation  for  Pt.  Here  S  denotes  the  sea  surface,  R  the  position  of  an 
area  dS  relative  to  the  origin,  and  n  the  normal  to  the  surface,  taken  as  positive  toward 
the  source.  We  will  take  the  origin  to  be  the  plane-surface  specular  point. 
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2.1  INTEGRAL  EQUATION  FORMULATION 


Since  5  is  a  pressure-release  surface,  one  boundary  condition  is  that  pj\s  =  0-  Then 
the  Helmholtz  integral  formula  (equation  5)  becomes 

Pt  (r)  =  po  (r)  (r)R)  ^ pr  (R) dS  .  (6) 

Two  integral  equations  [1]  can  be  developed  for  By  letting  r  approach  the  surface, 
we  have 


Po 


W  =  /  JsG{T\R)-g-pr(n)dS . 


(7) 


This  is  a  linear  integral  equation  of  the  first  kind  for  ^ .  Alternatively,  we  can  apply  the 
operator  =  n  •  V  to  equation  6  and  then  take  the  limit  as  r  approaches  the  surface. 
As  discussed  by  Meechan  [2],  it  is  legitimate  to  differentiate  through  the  integral  sign,  but 
the  limit  process  introduces  integrable  singularities  in  the  integrand.  The  result  is 


dpT 

dn 


dv0 
o  *  w 

~  dn 


-2JJsi  °wr>£»7'r) 


e 


dS  , 


(8) 


which  is  a  linear  integral  equation  of  the  second  kind  for  After  finding  ^  from 
either  equation  7  or  equation  8,  the  scattered  field  P  —  Pt  —  Po  follows  from  equation  6. 


Thorsos  has  compared  the  “exact”  results  of  numerical  solution  of  the  integral  equa¬ 
tion  to  solutions  obtained  by  using  the  Kirchoff  approximation  [3]  and  the  perturbation 
approximation  [4].  The  perturbation  approach  is  valid  when  the  root-mean-square  (rms) 
surface  height  h  is  small  compared  to  the  acoustic  wavelength  A.  Hence  it  is  applicable 
only  for  low-frequency  scattering. 


2.2  FRESNEL  CORRECTED  KIRCHOFF  APPROXIMATION 


To  obtain  the  Kirchoff  approximation,  begin  with  equation  8  and  drop  the  second  term 
on  the  right-hand  side 

dpr  _  ^dpo 
dn  dn 

Note  that  this  is  exact  for  a  flat  surface.  Hence  it  is  also  called  the  “tangent  plane” 
approximation.  Intuitively,  the  Kirchoff  approximation  would  appear  valid  when  the  sur¬ 
face  radii  of  curvature  are  large  compared  to  A  [5].  However,  Thorsos  has  found  the 
surface  correlation  length  /  to  be  the  most  important  parameter  in  determining  its  validity 
[3].  When  {  >  1,  the  Kirchoff  approximation  is  accurate,  except  at  low  grazing  angles 
where  shadowing  and  multiple  scattering  effects  enter.  Equation  6  may  therefore  be 
written  as 


o(rl 


f  [  —l „n  fP) 

J  J  dn  (‘  ”  ‘  ' 


e«qr-It| 


4x  |r  —  Il| ) 


R\  ^ 


no) 

v  ■*  -V 
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This  may  be  verified  by  first  noting  that  equation  10  gives 


'  *  //Hr+tH" 

=  //H^+§iG  dS' 


upon  use  of  the  boundary  conditions.  Then 


=  //[<*-'>§+(£+£) 
=  I  fa-P) 


G\  dS 


Finally 


/ /P"(R)|;G(r  |R,)<fS=jp(r) 
/  JC( r|R)  ^pb(H)<K--5P(r) 


The  last  expression  is  equivalent  to  equation  6. 

To  good  approximation,  the  direct  pressure  seen  at  a  point  r  resulting  from  a  spherical 
wave  source  located  at  rs  is  given  by 


Po  (r)  =  V-. 


2«fc|r— rs| 

lr  ~  rsl 


where  V  depends  only  on  direction.  As  discussed  by  Eckart  [5],  it  is  necessary  to  assume 
that  the  source  is  directional  and  ensonifies  only  a  finite  area  of  the  sea  surface.  The 
geometry  used  is  illustrated  in  figure  1. 

Now  if  |rs  —  R|  »  A  and  |r-  R|  >  A,  equation  10  may  be  further  simplified  to 


,(r)  =  // 


4w  |rs  -  R|  jr  -  R|  dn 


|e,fc(lrS-Rl+lr~Rl)|  dS  . 
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Figure  1.  The  scattering  geometry  used. 

For  a  sufficiently  directional  source,  the  denominator  may  be  replaced  by  its  value  at 
the  origin.  In  the  exponential,  we  will  use  the  Fresnel  approximation. 


|rs  -  R|  +  |r  -  R| 

\  =  rs  +  r  +  74'  (i,  y)  +  Qi  +  f3y  +  - bx 2  +  uxy  -f-  ^ej/2 

(19) 

where  C  is  the  local  surface  displacement,  a  zero-mean  random  variable,  and 

7 

=  —  (sin  V’s  +  sin  ip) 

(20) 

a 

=  —  (cost’s  —  COS  Xp  COS  ft) 

(21) 

0 

=  cos  4'  sin  ft 

(22) 

6 

=  —  sin2  rps  +  -  (l  —  c°s2  ^  cos2  ft) 

(23) 

r$  r  '  7 

e 

= - (-  -  fl  —  cos2  xp  sin2  ft) 

(24) 

r$  r  v  ’ 

V 

=  -  cos2  ip  cos  ft  sin  /i  . 

r 

(25) 
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The  importance  of  Fresnel  corrections  has  been  discussed  by  Melton  and  Horton  [6].  Note 
that 


dS  = 


dxdy  dxdy 


nz  cos  L 

where  V  is  the  angle  between  the  mean  normal  to  the  surface  (the  z-axis)  and  the  local 
normal  at  the  specular  point 

.  o2  +  /?2 


tan  T  = 


Thus 


jjj.  eiHrs+r) 

P  r)  =  - - 

rsr 


47 r  rsr  '  J  \  7  nz  7  nt ) 

exp  jzL  ( x,y )  +  ax  +  0y  +  ^ Sx 2  +  uxy  -f  ^ey2j  |  dxdy  . 


If  V  is  taken  as  constant  over  the  ensonified  area  and  a  Gaussian  surface-height  distri¬ 
bution  is  assumed,  this  can  be  evaluated  analytically  [7].  However,  a  physical  interpreta¬ 
tion  may  be  realized  if  the  integral  in  equation  28  is  evaluated  asymptotically  by  the 
method  of  stationary  phase  [8]. 


2.3  SPECULAR  POINT  THEORY 

Let  the  stationary  phase  points  be  denoted  Rj  =  (xj.y^Cj)-  A  second-order  expansion 
in  the  neighborhood  of  one  of  these  points  gives 

|rs  —  R|  4-  |r  -  Rj  =  rs  +  r  +  ax  +  f3y  +  ^<5x2  +  i/xy  +  ^ey2  + 

7  [Cj  L  (x  —  Xj)  dxCj  +  (y  —  yj)  ^yCjT- 

l  f29') 

-  (x  -  Xj)2  d2xC;  +  (x  -  Xj)  (y  -  yj)  d;  Cj  + 


\  (y  ~  Vi?  dl» 0 


That  the  phase  is,  in  fact,  stationary  requires 


0  =  7^iC;  +  a  +  6ij  +  vy$ 
0  =  idyQ  +  P  +  eyj  +  vxj 


5 


(32) 


a 


„  a  f>x}  +  t'V: 

:v>  ~ 


-  -?♦«(*) 
a,  o  =  -fi-'JU+oz. 

7  7 


=  -^  +  0 
7 


a- 


(33) 


Now  since  d*(  =  —  and  3VC  =  — jjf ,  it  can  be  seen  that  these  are  the  conditions  for 
specular  reflection,  a  result  that  might  have  been  anticipated.  Consequently 

|rs  -  R|  +  |r  -  Rj  =  rs  +  r  +  iQ  +  ax,  +  fiyj  +  +  vxjyj  +  ^cy]  + 

^7  [(*  ~  Xjf  dlx(j  +  2  (x  -  xj)  (y  -  Vj )  dxyCj+ 

+  j«  (*-*,')’  + 

v(x-Xj)(y-yj)  +  ^e(y-yjf  . 

Therefore,  p(r)  cm  be  written  as  a  sum  of  specular  point  contributions 


p(r)  = 


ik  eik{rs+r) 
4tt  rsr 


E7  1  + 


oM£' 

72  , 


Ij^i  , 


where 


and 


<t>j  =  k  (76  +  qxj  +  /?j/j  +  i^<5x*  +  uxjy}  +  \cyf) 
h  =  J  J  V exp  |  ~ik 7  [Aj  (x  -  x,)2  +  2R;  (x  -  Xj)  (y  -  y:)  + 

C,{y  -yj)2\}dxdy 


with 


a  -  #LCj  + 


=  «  +  - 
^  =  ^0  +  1. 


(34) 


(35) 


(36) 


(37) 


(38) 

(39) 

(40) 


If  V  is  taken  to  be  constant,  the  integral  can  be  evaluated  (8].  Three  different  situations 
must  be  considered,  corresponding  to  the  nature  of  the  phase  extremum  at  Rj .  For 
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convenience,  let  the  discriminant  Dj  =  AjC,  -  B]  ■  If  Dj  >  0  and  A,  <  0 ,  the  quadratic 
form  is  negative  definite  and  the  phase  function  has  a  relative  maximum.  If  Dj  >  0  and 
Aj>  0  ,  the  quadratic  form  is  positive  definite  and  the  phase  function  has  a  relative  mini¬ 
mum.  Finally,  if  D,  <  0 ,  the  quadratic  form  is  indefinite  and  the  phase  function  has  a 
saddle  point.  Then 


h  =  -Vti 


where  =  1,-1,  or  i  at  a  relative  maximum,  relative  minimum,  or  saddle  point,  respec¬ 
tively.  Thus 

p(r)  =  V - £  +  (42) 

2  VST  y  V  72  J  \Dj\t  1 

The  discriminant  can  be  related  to  the  Gaussian  curvature  K  of  the  surface  at  the 
specular  point.  Now 


K  = 


RaRb 


=  -  (€/)’]  [i  +  <a.o! + o-oT’  • 


where  Ra  and/4  are  the  principal  radii  [9].  Thus 


D,  = 


(RlJh)  j1  +  ^  (°2  +  P)  +  “  +  M  XJ  +  (ai/  +  W  + 

^  [(<52  +  u2)  x2  +  2 v  (6  +  e)  Xjyj  +  (e2  +  v2)  y]\  |  + 

\  +  «&<*  ~  2^0)  +  ~2  (&  -  v2) 

1  /  q2  +  02\2  „/Rj\ 

(RaRb),  [  +  72  )  +0(r)  • 


a2  4-  07 


Now  if  n  and  r2  see  the  same  specular  points,  a  good  approximation  if  they  are 
close  together,  we  have 

/  V  \  2  P>fc(ri-r2)  - -  ,, 

i  t  »  /  \  /  '  I  ^  V  ■»  *  //  r»  \  /  r»  —  /A.  \ 


/  V  \  eXK\r'-r*>  i -  , 

p(ri)p-(r2)  =  (— ) 


\2r5/  nrj 

/  ■p  N  2  etfc(r,-ra) 

\2r5/  rir2 


Z(R*R>h  + 


j#* 


At  this  level,  Fresnel  corrections  are  ignored. 


2.4  ENSEMBLE  AVERAGE 


The  phys’cal  interpretation  can  be  further  developed  if  the  ensemble  average  is  consid¬ 
ered 


(p(ri)p'  (r2)) 


ei/c(ri-r2) 


nr2 


X 


(£  (/J.  *)*<*■'<♦'-*■>> 


(48) 


Since  the  surface  is  random,  this  can  be  simplified  by  assuming  that  all  physical  quanti¬ 
ties  associated  with  the  specular  points  have  the  same  average  values  for  all  members  of 
the  ensemble,  although  the  number  and  location  of  the  points  fluctuate  from  member  to 
member  [8).  Therefore 

<£ (RM,)  =  +  =  <JV)< R.R„)  ,  (49) 


where  (N)  is  the  average  of  specular  points. 

In  the  analysis  of  the  second  term,  it  may  be  assumed  that  the  location  and  curvature 
of  a  specular  point  are  independent  random  variables.  This  allows  us  to  separate  the 
expected  values  of  the  phase  and  amplitude  factors.  Then,  if  the  phase  differences 
between  rays  from  different  specular  points  take  on  all  values  between  ±^m» r  with 
uniform  probability  [10],  the  desired  average  is 


<£  vispji Jh^UuhT^-^) 


(N)2(^RJib)2  -  (N)(RaRb) 


x 


(Sin  <f>max 

fimax 


(50) 


In  the  high  frequency  limit,  <j>max  »  1  and  this  term  may  be  neglected.  Thus 

/  V  \  2  pik^l-r-2) 

(P(ri)p-(rj))  =  If-)  - .  (51) 

\2rs /  r\r2 

We  will  briefly  mention  a  different  point  of  view1.  This  assumes  that  the  specular 
points  are  close  together,  so  that  the  phase  differences  are  approximately  equal.  Then 

(E  =  V  f  {N)2{y/KRby  -  (N)(RaRb)  ]  ,  (52) 

J*k  1  J 

for  some  constant  v  =  0(1),  and 

f  V  \2  e,<:(ri-rj) 

(p(r,)p-(r2))  =  (— ) 

'  Luby,  J.C..  and  D.  Funk,  private  communication.  The  Applied  Physics  Laboratory.  University  of  Washington, 

Seattle,  WA. 


(1  -  rj)  (N)(RaRb)  +  rj(N)2(^Rb): 


(53) 
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Note  that  the  dependence  on  geometry  is  implicit  in  (N),  since  (RaRb)  is  probably  not 
a  very  sensitive  function  of  the  geometry.  ( N )  will  be  a  maximum  near-normal  incidence, 
tps  =  §  .  As  *i>s  decreases,  more  points  fall  into  shadow  so  that  ( N )  decreases  and  will  be 
nearly  zero  at  grazing. 


2.5  PLANE-WAVE  SOURCE 

For  a  source  of  pb.ne  waves, 


Po  =  Pe,kS'r  . 


The  problem  is  that  of  scattering  from  ks  to  k .  Then 

P  d 


l(r'k  )  =  // 


|ci<:(kS'R+ir-Ri)  |  dS  . 


(54) 


(55) 


An  |r  —  R|  dn 

A  somewhat  different  expression  is  usually  encountered  in  the  literature  [8].  Observe  that 

r  =  [jr  —  R|2  +  P2  -f  2  (r  —  R)  •  R]  ^  (56) 

=  |r-R|  +  fc.R  +  0(^)  .  (57) 

Therefore,  ignoring  Fresnel  corrections, 

o'kr  r  r  Q 


p(r'k)  -  £Hpo;Mk(h'i)n}dS 

=  j  I P  — (kl-.k)  tHH-^dxdy  . 


The  remaining  analysis  proceeds  as  before.  Let 

(ks  -  k)  •  R  =  QX  +  /3y  +  7C  ( x,y ) 
and  expand  C  about  the  stationary  phase  points  to  get 

p(r.k)  =  ■ 


Finally 


1  „e»'*(ri-rj) 

(p(ri)p-(ra))  =  -P2 - (N)(RaRb)  . 


4  rjr2 

Note  that  the  scattering  cross  section  from  ks  to  k  is  therefore 

<r(ks,k)  =  Anr2  jp (r)|2  ~ 

=  n{N)(RaRb)  . 


(58) 

(59) 

(60) 

(61) 

(6?) 


(6?) 

(64) 


9 


Thus,  in  the  high-frequency  limit,  the  scattering  problem  has  a  simple  physical  interpreta¬ 
tion.  Each  specular  point  has  a  cross  section  equal  to  the  projected  aiea  of  a  sphere  whose 
radius  is  the  geometric  mean  of  the  principal  radii  at  the  point. 

Here  (N)  is  quite  sensitive  to  the  scattering  geometry,  independent  of  any  shadowing 
effects.  This  will  be  clarified  in  the  next  section. 


2.6  SURFACE  STATISTICS 


In  considering  the  effects  of  surface  statistics,  it  is  convenient  to  begin  with  a  surface 
roughness  spectral  density  W  (Kx,Ky)  normalized  such  that 


/oo  roo 

/  W  (A'r,  I\y)  dKxdKy  =  h2  , 
-oo  J  —  oo 


(65) 


with  h 2  =  (C2)  the  mean  square  surface  height  and  I\x  and  Ky  the  surface  spatial 
wavenumbers.  The  roughness  spectral  density  is  related  C  to  by 

|2 


W(I<X,KU)  =  lim 


2  Wo 


x.-*~  (2n)  I2 


Ft  f\L  a*,y)e-i{I<*x+K'>v)dxdy 
J-hL J-\L 


) 


(66) 


2.6.1  Gaussian  Roughness  Spectrum 

Typically,  a  Gaussian  roughness  spectrum  is  used 


An 


W  (I\x,  Ky)  =  —  e-H^+A»)': 


(67) 


where  /  is  again  the  surface  correlation  length.  The  corresponding  correlation  function  is 

/OO  too 

/  W  (A'x,  A'y)  e'^x+h^dKxdI<y  (68) 

-OO  J  —  oo 


-oo  J  —oo 


(69) 


Barrick  [11|  has  derived  expressions  for  the  average  number  of  specular  points  per 
unit  area,  ua,  and  the  average  curvature  in  terms  of  the  surface  statistics,  for  plane-wave 
illumination.  For  an  isotropic  surface  with  Gaussian  statistics 


7.255  _u 4x 


(70) 


where  s2  is  the  mean  square  value  of  the  total  slope  at  a  point  on  the  surface, 
»2  =  <(C«)!  +  (0,<)’>  =  4  +  »;•  Here 


,s2  =  4 


l 2  ' 


(71) 
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r  again  is  the  angle  between  the  mean  normal  to  the  surface  and  the  local  normal  at  the 
specular  point.  The  average  curvature  can  be  approximated  as 


(K)  = 


7.255  ,s2 
7 r  / 2 


cos4  r 


so  that 


{RaRb)  =  0.1378tt 


/2 


sec4  T  . 


Therefore,  the  scattering  cross  section  per  unit  area  is 


(72) 


(73) 


(74) 


2.6.2  Pierson-Moskowitz  Roughness  Spectrum 


Fully  developed  sea  surfaces  exhibit  a  range  of  roughness  scales  (multiscale  surfaces). 
However,  the  Gaussian  spectrum  describes  a  single-scale  surface.  A  more  realistic  model 
is  that  of  a  Pierson-Moskowitz  surface,  in  which  the  spectrum  is  completely  determined  by 
the  wind  speed.  To  obtain  the  spatial  spectrum,  begin  with  the  Pierson-Moskowitz  fre¬ 
quency  spectrum  [12]  given  by 


with  a  =  S.10  x  10  3,  /?  =  0.74,  g  =  9.81  m/s~2  and  U  the  wind  speed  at  a  height  of  19.5 
m.  The  gravity  wave  dispersion  relation, 

u>J=*|A'|,  (76) 

may  be  used  to  relate  .5(w)  to  a  “nondirectional”  wavenumber  spectral  density  [13]. 


W(A')  =  SM  £.=S (vSa)  jv® 


ir  K  >0,  and  VV (I\)  =  0  for  A'  <  0 .  This  is  normalized  such  that 

fwW  =  h 2  . 

Jo 


(77) 


(78) 


Thus 


W(A')  = 


(79) 
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and 


=  \f-^-u2 

V  w 

=  (5.32  x  IQ'3  m)  . 


(80) 

(81) 


The  definition  of  the  2D  wavenumber  spectrum  W(Kx,Ky)  requires,  in  addition 
to  equation  75,  a  model  for  the  azimuthal  dependence.  Scattering-theory  con¬ 
ventions  require  that  W  (K)  =  VT(-K).  A  transformation  to  polar  coordinates,  using 
W  (Kx,  I\y)  dKxdI<y  =  Wi  (A',  <t>)  KdKd<f> ,  gives 


W  (Kx,Ky)  =  Wx  (jKl  +  Kl,  arctan  h±\ 


(82) 


If  the  wind  direction  is  chosen  to  be  in  the  x  direction,  then  4>  is  also  the  azimuthal  angle 
relative  to  the  wind.  Then 


W  ( K  ^ 

=  — jr-l,  (83) 

where  describes  the  azimuthal  dependence.  It  is  normalized  such  that 

/%(A',<£)<ty=l  ,  (84) 

Jo 

and  satisfies  the  symmetry  property  $  ( A",  <£  +  7r)  =  $  (I\,(f>).  A  cos2  azimuthal  depend¬ 
ence  for  W  is  often  assumed,  so  that 

$  (k,  </>)=  —  cos2  <j> .  (85) 

7T 

In  this  case 


W(I<x,I<y) 


Q  Kl  (  Pg2  1  \ 

(a'|  +  k;)3  exp  V  + 


(86^ 
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3.0  IMPLEMENTATION  CONSIDERATIONS 


3.1  POINT-SCATTERER  MODEL 

For  the  purpose  of  simulation,  the  point-scatterer  model  appears  the  best  implementa¬ 
tion  of  forward  scattering.  This  is  because  specular  point  theory  gives  reasonable  fidelity 
without  excessive  computational  loads.  Autoregressive  and  inverse  beamformer  models, 
although  computationally  efficient,  are  poor  representations  of  sources  that  are  narrow- 
band  or  have  significant  Doppler  spread.2  They  are  also  less  easily  visualized. 

The  model  is  also  suitable  for  the  analysis  of  short  pulses.  The  times  at  which  specular 
points  “turn  on”  and  “turn  off”  can  be  calculated,  and  consequently  the  amount  of  time 
required  to  ensonify  the  region  required  for  steady-state  analysis.  For  short-enough  pulse 
lengths,  the  steady-state  condition  is  not  reached.  Therefore,  a  higher  coherence  would  be 
observed  than  that  predicted  by  steady-state  analysis  [11]. 

The  most  straightforward  implementation  begins  with  the  partition  of  a  sufficiently 
large  portion  of  the  surface  into  a  grid  with  spacing  of  a  correlation  length.  For  each 
correlation  area,  nA  is  used  to  determine  whether  a  specular  point  is  present.  If  so,  it  is 
randomly  placed  in  the  area.  These  determine  the  paths,  which  are  then  summed.  Luby 
and  Funk  [11]  have  compared  this  technique  with  the  closed-form  solution  of  McDaniel 
[7],  finding  good  agreement  except  when  the  source  and  receiver  depths  differ  signifi¬ 
cantly. 


3.2  RANDOM  SURFACE  GENERATION 


Single-ping  realizations  may  be  analyzed  by  generating  random  surfaces  for  specified 
roughness  parameters  and  then  determining  the  specular  reflections.  There  are  several 
good  reasons  for  doing  so.  Examination  of  ping-to-ping  fluctuations  provides  a  measure 
of  higher-order  statistics  unavailable  from  ensemble  methods.  These  are  expected  to  be 
significant  in  the  performance  of  adaptive  beamformers.  In  addition,  the  actual  simulation 
of  surfaces  provides  a  tool  for  the  validation  of  methods  proposed  to  calculate  coherence. 


Consider  a  surface  that  is  random  in  both  directions.  The  objective  is  to  find  ((x,y) 
at  NxM  points  with  spacing  A.t  and  Ay  over  lengths  LT  =  NAx  and  Ly  =  MAy. 
Realizations  can  be  generated  at  points  xn  ~  nAx  and  ym  =  mAy  by  using 


C  (^ni  J/m  )  — 


LXL 


i 


2  1 


(37) 


2  Hidinger,  R.  M.,  and  J.  M.  Zyphur,  private  communication.  Naval  Ocean  Systems  Center,  San  Diego,  CA. 
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where 


F(Kj,  I<i)  =  2ir  [LxLvW  (\Kj\, |A',|)J*  G}G,  . 


For  j  >  0  and  l  >  0, 


Gj  = 


Gi  = 


( 

( 


N(  0,1) 

N(  0,1) 

J-[iV(  0,1)  + .TV  (0,1)] 


if  j  =  0  or  y 
otherwise 

if  /  =  0  or  f 
otherwise 


(88) 


(89) 

(90) 


while  for  j  <  0,  Gj  =  Gi, ,  and  similarly  for  /  <0,  Gi  =  Gl, .  Here  K,  =  xf,  Ki  —  , 

and  each  occurrence  of  N  ( 0,1)  indicates  an  independent  sample  taken  from  a  zero-mean, 
unit-variance  Gaussian  distribution. 

For  numerical  implementation,  it  is  convenient  to  rewrite  equation  87.  First  note  that 


Um  — 


1 


1 


£-1 
2  1 


f-i 


AxA y  N M 


E  £  FiK^Ki)  ***»*') 


(91) 


M 

2 


Shifting  indices  by  letting  *7  =  j  +  y  and  £  =  /  +  y  gives 


-  — 7r;(n+m)  i  N-lAf-1  v  /  x 

e  .  1  y;  rf  /{  n,a'c  „)e2!"(*J+*£) 
AxAy  NM  v 

(-l)n+m 


Ax  Ay 


-ifft2  (F)  , 


(92) 

(93) 


where  ifft2  denotes  the  two-dimensional  inverse  Fourier  transform.  This  is  an  efficient 
form  for  calculation. 

Figure  2  shows  an  example  of  a  surface  realization.  Here  the  grid  was  taken  to  be  128 
x  64  with  spacings  Ax  =  Ay  =  0.2  m,  The  wind  speed  was  chosen  as  U  =  10  m/s-1.  A 
slice  taken  along  the  x  direction  of  this  surface  is  shown  in  figure  3.  Some  of  the  small- 
scale  variation  results  from  the  relatively  coarse  resolution  used.  The  latter  figure,  in 
particular,  indicates  the  multiscale  behavior  of  the  surface. 
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zeta  (m) 


5 


4.0  SUMMARY 


Forward  scattering  should  be  best  simulated  with  a  point-scatterer  model  through  the 
use  of  specular  point  theory.  A  derivation  of  specular  point  theory  from  the  Helmholtz 
integral  formulation  was  presented,  and  the  attendant  approximations  noted.  All  hold 
reasonably  well  for  most  cases  of  interest,  and  produce  only  minor  errors  with  respect  to 
the  “exact”  results. 

Several  implementation  issues  were  discussed  as  well.  A  simple  way  to  realize  the 
point-scatterer  model  is  through  the  partition  of  the  surface  into  correlation  areas  and  the 
probabilistic  assignment  of  specular  points.  Detailed  examination  of  ping-to-ping  fluctua¬ 
tions,  or  the  effects  of  higher-order  statistics,  require  the  generation  of  surfaces  with  a 
desired  roughness  spectrum.  A  procedure  to  realize  this  was  presented. 
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